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We report on the structure and adsorption properties of Pd^ (n = 1 — 4) clusters 
supported on the rutile Ti02 (110) surfaces with the possible presence of a surface 
oxygen vacancy or a subsurface Ti-interstitial atom. As predicted by the density 
functional theory, small Pd clusters prefer to bind to the stoichiometric titania surface 
or at sites near subsurface Ti-interstitial atoms. The adsorption of Pd clusters changes 
the electronic structure of the underlying surface. For the surface with an oxygen 
vacancy, the charge localization and ferromagnetic spin states are found to be largely 
attenuated owing to the adsorption of Pd clusters. The potential energy surfaces 
of the Pd monomer on different types of surfaces are also reported. The process of 
sintering is then simulated via the Metropolis Monte Carlo method. The presence of 
oxygen vacancy likely leads to the dissociation of Pd clusters. On the stoichiometric 
surface or surface with Ti-interstitial atom, the Pd monomers tend to sinter into 
larger clusters, whereas the Pd dimer, trimer and tetramer appear to be relatively 
stable below 600 K. This result agrees with the standard sintering model of transition 
metal clusters and experimental observations. 
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I. INTRODUCTION 



In the recent decades, the heterogeneous catalysis of metal/metal-oxides has advanced to 
the sub-nano regime, where one promising catalyst architecture is metal clusters composed 
of only a few atoms and supported by oxide surfaces. The chemical properties of these sub- 
nano metal clusters have been shown to be very different from those of corresponding bulks 
or nanoparticles with hundreds of atoms. In model studies of chemical reactions catalyzed by 
these systems, the oxides surface does not only provide the support for the cluster but also 
directly participates in the reaction by supplying charges, adsorption sites, and additional 
reaction paths. As one of the most popular model systems on the transition metal-oxide 
surfaces, the rutile Ti02 (110) provides rich and controllable surface conditions for studies 
on various catalytic reactions.'^' 

Experimentally, though the Ti02 surface is often prepared in ultra high vacuum (UHV) 
conditions to maintain the cleanliness of the surface as much as possible, it is known that 
many Ti02 surfaces prepared by ion sputtering and annealing in vacuum have localized 
band-gap states with Ti3rf character.!^ These gap states, which lie about 0.8-1.1 eV below 
the conduction band minimum (CBM), are believed to be involved with imperfections of the 
surface, such as surface oxygen vacancieJ^^^ or near-surface Ti-interstitial atoms.'^^'^ What 
roles these defects play in the metal/metal-oxides catalysis and how they interact with the 
metal clusters is not well-understood. 

On the scale of several atoms, the size of the metal clusters can have a dramatic effect 
on the catalytic activity.'^^^ For example, in a recent temperature-programmed reaction 
experiment done by the Anderson group, the carbon monoxide oxidation rate in vacuo 
shows a 150% jump when the size of the deposited Pd^ cluster increases from monomer 
(Pdi) to dimer (Pd2).^^ In the meantime, it is well-known that metal clusters are vulnerable 
to sintering processes, in which small sized clusters gradually group into larger ones and 
finally lose their unique catalytic activity.'^SHZS 

A good understanding of these intriguing catalytic behaviors of Pd clusters on titania 
requires an extensive microscopic knowledge of how Pd clusters bind and move on the 
titania surfaces with the possible existence of surface defects. A theoretical study, with 
controlled surface cleanliness, particle size, and number density, is ideal for this purpose. In 
this paper, based on the density functional theory (DFT) with Hubbard U correction, we 
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present a thorough study on the adsorption geometries, energetics, electronic structures, and 
magnetic properties of Pd^ (n = 1 — 4) clusters as well as the potential energy surface (PES) 
of the Pd monomer on the stoichiometric and defective titania surfaces. With these data we 
perform Monte Carlo (MC) simulations and further investigate the sintering patterns of Pd 
clusters under experimentally relevant conditions. 

II. METHOD OF CALCULATION 

The calculations were performed at the DFT level with a plane-wave basis set and ultra- 
soft pseudopotentials implemented in Quantum Espresso?^ The Ti 3s- and 3p-states were 
treated as valence states. The number of valence electrons considered for Ti, Pd, O and H 
were 12, 10, 6 and 1 respectively. The cutoff energy for the plane- wave basis set was chosen 
to be 32 Ry 435 eV). When the surface oxygen vacancy or Ti-interstitials are present, the 
additional electrons from these defects may localize on nearby Ti atoms and, reduce these 
Ti^+ ions to Ti^+. These highly localized states are known to be poorly described by tra- 
ditional density functional, including the local density approximation and the generalized 
gradient approximation (GGA). To describe the exchange-correlation effect of valence elec- 
trons, we adopted the DFT+[/ approach (implemented in the simplified DudareV^ scheme) 
with the Perdew-Burke-Ernzerhof (PBE) functional.^ The U term was only applied to Ti 
atoms with its value being computed self-consistently^'' as 4.2 eV. Spin-polarized calcula- 
tions with fixed multiplicity were performed in all cases. A gaussian smearing with width 
of 5 mRy 68 meV) was applied to the occupation of electron bands. The space charges 
were assigned to individual atoms through the Bader charge analysis .'^^'^ 

With above settings, the lattice constant of titania in equilibrium reads a = 4.67 A and 
c = 3.02 A, which agrees with the experimental value (a = 4.58 A and c = 2.95 A) with 
errors smaller than 2%. These values were used to construct a stoichiometric (4 x 2) surface 
slab, which was the starting point of all subsequent calculations with adsorbates. One Pd 
cluster per supercell gives 1/8 ML coverage, where 1 ML is equivalent to the density of 1 
cluster per Ti02 (1 x 1) unit cell (5.2 x 10^^ cm~^). The vacuum distance between two 
images of slabs is chosen to be 13 A, contributing an error in energy of the order of 2 meV. 
The size of the supercell is large enough to use a single point (shifted from the F point) in 
the reciprocal space to evaluate the total energy. The adsorbates were only placed at one 
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side of the slab. The energies of gas phase Pd^ cluster were calculated in a 15 A cubic box at 
the r point. A good balance between the computation speed and accuracy can be achieved 
using a four-trilayer of Ti02 unit with the top two trilayers fully relaxed and the remaining 
two constrained to their bulk equilibrium positions. To avoid the band-gap narrowing and 
spurious surface state owing to this additional constraint, we saturate the dangling bonds 
of Ti and O at the back of the slab using pseudo-hydrogen atoms with fractional charges 
4/3 and 2/3 respectively, as suggested by Kowalski et all^^^ The total number of atoms 
including the pseudo- hydrogens per stoichiometric supercell is 208. 



III. RESULTS AND DISCUSSION 



In this section, we show results on the adsorption properties of Pd^ clusters (n = 1 — 4) 
on the stoichiometric and two defective Ti02 (110) surfaces: the surface with an oxygen 
vacancy and the one with a subsurface Ti-interstitial atom. With the energetics obtained 
from DFT+U calculations, we perform MC simulations on the growth and sintering of small 
Pd clusters. 

To facilitate the characterization of major adsorption sites, we introduce the following 
notions (taking the stoichiometric surface as an example), shown in Fig.[l] The five- and six- 
fold coordinated Ti atoms are marked as Ti^f and Tie/. The two-fold coordinated surface 
oxygen atom in the bridging row is labeled as O5, whereas O5 represents the three-fold 
oxygen atom in the first trilayer of the Ti02 unit. The hollow site between the two O5 and 
a nearby O5 atom is denoted as h^. 



A. Gas phase Pd^ (n = 2 — 4) cluster 

The Pd atom has a M^^bs^ configuration. The sp-d hybridization occurs when Pd atoms 
are brought into contact, leading to partial occupancy in the s- and d- atomic orbital. The 
exchange energy among these partially filled 4rf states tends to drive the system to higher 
multiplicity. In bulk, though being paramagnetic, Pd is very close to satisfy the Stoner 
criterion^^ and, becomes ferromagnetic with a small lattice expansion. I^^El 

Several studies have been dedicated to the energetics, structures and magnetic properties 
of Pd clusters in gas phase.™! In the case of small Pd clusters (n = 2 — 4), the ferromagnetic 
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FIG. 1. The stoichiometric mtile Ti02 (110) surface. Atoms are labeled as follows: Ti^f and Tigj 
— five- and six-fold coordinated Ti atom; O5 — surface bridging oxygen atom; O5 — first layer 
in-plane oxygen atom; h^, — hollow site between two O5 atoms and a nearby O5 atom. 

configuration with magnetic moment M = 2 /x^ was identified to be the global energy 
minimum. The spin-singlet isomers are generally a few tenths of eV higher in energy. In 
Tab. [l| we give the binding energy per atom: E}j{n) = — [£'(Pd^) — n£'(Pdi)]/n, where 
E(X.) is the total energy of the component X, the dissociation energy: Ed{n) = £'(Pd^_i) + 
£'(Pdi) — £'(Pd^), which gives the energy cost of removing one Pd atom from the Pd^ 
cluster, and the average Pd-Pd bond length for n = 2 — 4. The binding energies obtained 
by Futschek and co-workerd^ are also listed for comparison. 

In general, the average bond length of Pd cluster at M = 2 /x^ is significantly smaller than 
the bulk value 2.8 A). As n increases, d and £"5 increase accordingly because the average 
number of dangling bonds per atom decreases. For the n = 4 case, though the tetrahedral 
structure has significantly lower energy (0.81 eV) in gas phase, the square structure becomes 
important when the underlying symmetry of the Ti02 surface is taken into account. The 
dissociation energy of Pda and Pd4 are much larger than that of Pd2, suggesting in gas phase 
the dimer is much easier to break into separated entities (monomers). 

For the low energy spin-triplet states, we found the highest occupied molecular orbital 
(HOMO) of Pd3 is degenerate and partially occupied. In this case, introducing the Jahn- 
Teller distortionP^ can break the symmetry, and further lower the energy. Hence these 
obtained triplet states are not the true ground states. As pointed out by Futschek et al^ 
the GGA functional is not ensemble v-representabld^ and hence incapable of describing this 
effect correct ly.l^ However, we argue that the introduction of the surface-cluster interac- 
tion can effectively break the symmetry and remove such degeneracy, leaving only the gas 
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n Geometry M (hb) d (A) (eV) Eb (eV) [Ref. 37J E^ (eV) 



2 Dimer 


2 


2.51 


0.634 


0.646 


0.634 







2.59 


0.505 


0.473 


0.505 


3 Triangle 


2 


2.55 


1.236 


1.250 


2.442 







2.52 


1.178 


1.250 


2.267 


4 Tetrahedron 


2 


2.63 


1.648 


1.675 


2.884 







2.63 


1.584 


1.654 


2.626 


4 Square 


2 


2.52 


1.445 


1.485 


2.071 







2.51 


1.416 


1.234 


1.954 


Bulk 





2.81 


3.671 


3.704 





TABLE 1. Geometries, magnetic moment M, average bond length d, binding energy £"5 and 
dissociation energy of the Pd^ clusters for n = 2 — 4. The corresponding value of Pd bulk is 
also given for comparison. 

phase results affected. This gives a constant shift in the calculated adsorption energies. No 
qualitative change of our main conclusions is expected. 

B. Pd^ on the stoichiometric surface 

The stoichiometric Ti02 is an insulator with a band-gap of 3 eV. Its (110) surface can 
be easily prepared in UHV conditions.'^ As shown in Fig. [!} the surface can be viewed as 
a composition of alternating rows of Tisj atoms and the bridging O5 atoms, both of which 
have dangling bonds owing to the surface cleavage. 

We compute the adsorption energy as the difference of DFT total energies of the two 
components (the cluster and the surface) before and after the adsorption: 

EaM = E[Smi\ + E[Fdn] - £;[Surf+Pd,], 

which is always positive by definition. The sintering energy Es[n] has a similar definition of 
the dissociation energy in gas phase (see Sec. |III A ): 

Es[n] = £;[Surf+Pd^_i] + £;[Surf+Pdi] - £;[Surf+Pd^] - E[Smi\. 

Eg measures how easily a Pd^ cluster can be decomposed into a Pd^_i and a Pd monomer 
on the surface. The adsorption geometries for Pd^ with the highest adsorption energies are 
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M=0; E,=\.5?>qW 




(b) 



2.70 



2.51 




M=0; £,^=2.33eV 



2.55 (d) 




M = 2\ =2.48eV 



M=0; £,^=2.41eV 



FIG. 2. Adsorption geometries of Pd^ cluster for n = 1 — 4 on the stoichiometric Ti02 (110) 
surface. The corresponding magnetic moment M and adsorption energies E^a are given below the 
geometries. Selective bond lengths (in A) of Pd-Pd and Pd-0 bond are given with arrows. Insets: 
ground state geometries of Pd^ cluster in gas phase. In the case of Pd4, the ground state is the 
tetrahedron whereas the square planar structures has energy 0.81 eV higher. 

shown in Fig. [2j The adsorption energies, sintering energies, magnetic moment and average 
transferred electrons of selected atoms (including Pd and the nearby O5 atom bonding to 



Pd) are listed in Tab. II 



The studies of the adsorption properties of Pd clusters on stoichiometric Ti02 surface^^^^SH 
are not as abundant as those of other transition metal clusters on the titania surfaces. The 
Pd monomer and dimer were investigated by Bredow and Pacchioni, who showed that the 
preferred binding site for Pdi is between two O5 atoms with calculated adsorption energy 
1.0 eV.'^ Sanz and Marquez argued that allowing relaxation of the surface gives the h^ as 
the preferred adsorption site of the Pd monomer, which is inclining towards a nearby O5 
atom and forms a right angle with Tisj below.'^^'^ Our results confirms Sanz and Marquez's 
conclusion (see Fig. [2] (a)). We found that adsorption of the Pd monomer on the h^ site is 
0.3 eV lower in energy than that between two O5 atoms. 

The Pd dimer has several configurations of adsorption in a small energy range. One 



n Geometry M {jj^b) Eadie^) £^s(eV) Ag[Pd](e/atom) AQ[05](e/atom) 



1 Monomer 





1.53 




-0.24 


0.07 


2 Dimer 





2.33 


0.54 


-0.16 


0.06 


3 Triangle 


2 


2.48 


1.06 


-0.13 


0.05 


4 Square 





2.41 


1.28 


-0.11 


0.06 


4 Tetrahedron 


2 


2.40 


1.27 


-0.08 


0.05 



TABLE 11. Magnetic moment, adsorption energy Ead^ sintering energy Eg^ and electron transfer 
AQ of Pd^ (n = 1 — 4) clusters adsorbed on stoichiometric Ti02 surface. The minus sign in AQ 
means a decrease in the number of electrons on the atom. 

possible way is to adsorb on two adjacent h^ sites along the O5 row, as discovered by 
Sanz and Marquez. In this case, the distance between the two Pd atoms in the dimer is 
stretched from 2.51 A to 2.75 A, to commensurate with the underlying Ti02 lattice (the 
lattice constant is 3.06 A along the (001) direction). Another configuration, which is 0.12 
eV lower in energy, is the Pd dimer bridging two rows of O5 atoms with an angle of about 
120 degrees, as shown in Fig. [2] (b). Both Pd atoms are still on the h^ sites. The sintering 
energy, Eg^ which is the energy cost of breaking the dimer into two monomers, is 0.54 eV. 
This value is slightly lower than that in the gas phase (0.63 eV), suggesting that the Pd 
dimer is relatively stable on the stoichiometric titania surface. 

For the Pd trimer and tetramer, one notices that because of the cluster-surface interaction, 
the adsorbed structures are all in-plane planar. However, for the Pd4 cluster in gas phase, 
we found that the tetrahedron Pd4 is 0.81 eV lower in energy than the square Pd4 isomer, 
whereas the two isomers have almost identical energies when adsorbed on titania. The square 
structure is slightly more favored (by 0.01 eV). The adsorption geometry of the square 
structure is shown in Fig. [2] (d). The incommensurability between the surface structure 
and the ground state cluster geometry in the gas phase forces the Pd-Pd bond length to be 
stretched by about 6% on average. The sintering energies of Pd3 and Pd4 are lower compared 
to their gas phase values owing to the cluster-surface interaction. The tendency of wetting 
the surface is quite different from the case of the noble metals, such as Au and Pt, the clusters 
of which prefer to grow along the surface normal, leading to either 3D structures, or out-of- 
plane "standing" geometries J ^^ ' ^^ ' The in-plane adsorption picture confirms the experimental 
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observation of the Anderson group, of which the authors conclude that for n < 10 the 
adsorbed Pd clusters are mostly in-plane 2D. The hollow site h^ is consistently the most 
preferred adsorption site for the Pd atoms, as long as the geometry of the cluster allows. 

In general, we found that the energy differences between the spin-singlet and the spin- 
triplet states are in the range of 0.2-0.3 eV, though their relative order may interchange 
compared to the gas phase results. For example, the cluster-surface interaction quenches 
the high-spin state of the Pd dimer, whereas the spin-triplet states of the Pd trimer survives. 
To characterize the oxidation state of the adsorbed Pd^ clusters, we performed Bader charge 
analysis to quantify the amount of charge transfer between the surface and the cluster, which 
is shown in Tab. [11} The average number of electrons on the surface O5 atom before the 
adsorption reads 7.04. We see that there is no significant charge transfer between the cluster 
and the surface. The bridging oxygen atom which forms bonds with Pd cluster gains 0.06e on 
average. Compared to the charges in gas phase, Pdi donates about 0.24e to the surface, while 
for Pd2-Pd4 these values increases to 0.32e, 0.39e and 0.44e respectively. The transferred 
charges are not evenly distributed within the cluster. In general, this anisotropic distribution 
of charges could infiuence the surface chemistry, in particular, the adsorption and activation 
of certain gas molecules. This aspect will be further explored in later publications. 

The small charge transfer indicates that Pd clusters are only slightly oxidized, largely 
retaining their gas phase chemistry. To further understand the nature of the surface-cluster 
interaction in terms of molecular bonding, we computed the density of states (DOS) of the 
supercell and plotted selected Pd states inside the Ti02 band-gap, shown in Fig. [3| For Pdi, 
the original degenerate 4d states now split into five states spanning an energy window about 
1.3 eV. Though these gap states largely retain the rf-character, orbital hybridizations with 
nearby O5 and O5 atoms are visible. For n = 1 — 4, all HOMO states suggest weak bonds are 
formed between the Pd atoms and the nearby surface O5 and O5 atoms. By performing the 
so-called constrained space orbital variation method,!^ Sanz and Marquez argued that the 
surface-cluster interaction comes largely from orbital polarization of Pd,'^ which can be seen 
from the spatial distribution of the HOMOs. The relative height difference in the DOS of 
the band-gap states comes from the symmetry requirement imposed by the underlying Ti02 
surface. For the n = 4 case, the adsorption of square-shaped cluster is slightly preferred over 
that of the tetrahedral one, because the planar cluster has a large number of states with 
four- fold symmetry imposed by the surface (shown as a high peak in the DOS), as the inset 
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E(eV) 

FIG. 3. Density of states (DOS) of the adsorbed Pd^ clusters (n = 1 — 4) on stoichiometric 
titania surface. The DOS of the pure Ti02 surface is plotted as black dashed line. For each n, 
selected gap-states are shown in the insets. The cut off of the isosurface is chosen to be 0.001 
e/bohr^. (The cut off is fixed throughout the paper unless stated otherwise.) The arrows mark 
their corresponding positions on the DOS. 

of Fig. [3] (d) shows. This effect makes the planar cluster having the energy comparable to 
that of the tetrahedral structure, completely closing the gap of 0.81 eV between the ground 
state energies of the two geometries in gas phase. 

C. Pd^ on the surface with an oxygen vacancy 

Surface vacancies are the most common point defects seen on Ti02 surface. When the 
bridging oxygen is removed, a pair of electrons is left to fill the nearest available states. 
Experimentally, the spectrum corresponding to Ti^+ states were observed by X-ray photo 
spectroscopy (XPS) when oxygen vacancies were created, suggesting the existence of local- 
ized TiSrf state from the reduction of the electron lone pair. Traditional DFT methodpl^ 
failed to capture this phenomenon owing to the well known self-interaction errors^^ built 
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in the approximate functionals. Using the hybrid B3LYP functional, Di Valentin et al 
successfully obtained the gap states for systems with oxygen vacancies. The DOS shows 
an energy split of about 0.1 eV in the Ti02 band-gap. Within the framework of DFT+[/, 
Morgan and Watson showed that there are two different regimes describing the localization 
of electrons. For U < 4.2 eV, electrons partially localize on several subsurface Ti ions. A 
qualitative change was seen when U was increased above 4.2 eV. 

With the self-consistently calculated U (4.2 eV) on all Ti atoms, we found that multiple 
sites of localization can be achieved, agreeing with Deskins et a/.'s conclusion.'^ The lowest 
energy configuration involves one electron localizing on one Tigj neighboring the vacancy 
and one on a next nearest neighbor subsurface Ti ion under the surface Ti^f row. The 
spin-triplet ground state has splitted DOS in the titania band-gap, as shown in Fig. |4} The 
localization picture is similar to Calzado and co-workers' results.'^ Owing to the inequivalent 
localization sites, both energies and orbitals of the two gap-states differ from each other. 
The electron localizing near the surface vacancy is lower in energy (by 0.6 eV). 

The reason that various localization sites were obtained roots deeply in the use of the 
parametric approach of correcting the self-interaction error of density functionals. Different 
choices on the values of U or the percentage of the Hatree-Fock exchange included in the 
hybrid functional determines not only the energy of the gap states, but also the positions of 
the localized electrons.^'' (see Ref. [58] by Dohnalek and co-workers for a complete review.) 
Experimentally, based on the resonant photoelectron diffraction measurement, Kruger et al 
found that the excessive charges are largely shared by several surface and subsurface Ti 
sites with a dominant contribution from the subsurface Ti ioiF^, which is in line with our 
calculation (see Fig. On the other hand, using the low temperature STM technique, 
Minato et al showed that the two excessive electrons do not localize directly on the vacancy 
but on neighboring surface Ti^f ions.^° In reality, the dynamics of the surface structures at 
non-zero temperatures may play an important role in the electronic localization. Combining 
DFT+U and the Car-Parrinello molecular dynamics, Kowalski and co-workers have shown 
that the localized electrons are in fact very mobile, and "hop" between different sites because 
of the thermal vibration of the Ti02 lattice.l^ Both splitted and non-splitted DOS were 
obtained depending on the transient state of the localized electrons. The subsurface layer 
was found to be occupied by excessive electrons about 70% of time, with the surface layer and 
the third trilayer populated at 20% and 10% of time respectively. This result suggests that 
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FIG. 4. DOS of the defective Ti02 surfaces with an oxygen vacancy. The inset shows the structure 
of the corresponding surface. The position of the locahzed electron is shown in golden color. The 
spin-resolved DOS are shown as solid black and dashed red lines for spin up and down respectively. 
The two gap states with localized electrons are shown in the insets with arrows marking their 
positions on the DOS. 



most of previous works (including this paper) actually captured snapshots of this complicated 
phenomenon at zero temperature. It would be very interesting to see how the adsorbates 
interact with the localized electrons near the oxygen vacancy at finite temperature. In the 
following we show that, even at zero temperature, the adsorption of additional molecules 
near the oxygen vacancy can affect the local geometry and hence change the energetics and 
electron localization sites in various ways. 

Fig. [5] (a-d) shows the adsorption geometries and energetics of a Pd^ cluster close to a 



surface oxygen vacancy. The magnetic moment and sintering energies are listed in Tab. Ill 



In Sec. |IIIB , we show that on the stoichiometric surface, Pd has a tendency to binding to 



protruding O5 atoms. This picture is further confirmed in the adsorption nearby a vacancy. 
Contrary to the results of Ref . UK the preferred adsorption site for Pdi is not on the vacancy 
but away from it. And the "usual" h^ site next to an atom now is less favored (by 0.25 
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n Geometry M (/is) Ead{(i\) Es{eY) AQ[Pd](e/atom) 



1 Monomer 


2 


1.79 




-0.28 


2 Dimer 





2.36 


0.04 


0.26 


3 Triangle 





1.99 


0.02 


0.13 


4 Tetrahedron 





1.37 


0.22 


0.10 



TABLE III. Magnetic moment, adsorption energy Ead^ sintering energy Eg^ and electron transfer 
per Pd atom AQ of Pd^ (n = 1 — 4) clusters adsorbed on the Ti02 surface next to an oxygen 
vacancy. 



eV) as compared to the bridging site between two O5 atoms (see Fig. [5] (a)). The total 
energy gain by moving away from the vacancy to the bridging position between two O5 
atoms is 0.65 eV. The adsorption energy of the geometry shown in Fig. [5] (a) is 0.26 eV 
higher than that of the Pd monomer on stoichiometric surface. The presence of the Pd 
monomer also affects the localization site of the two excessive surface electrons — both of 
which now localize in the subsurface layer, in an effort to maximize their mutual distance 
and those with the Pd monomer. Such local repulsion effect cannot be captured without 
correcting the self-interaction error in the density functional. 

Similar to the results in Sec. Ill B| the Pd dimer has multiple adsorption geometries 



within the energy range of 0.2 eV. Fig. [5| (b) shows the configuration with the lowest energy, 
the Pd dimer shows no preference over the type of surface — either stoichiometric or with 
an oxygen vacancy, the adsorption energies are roughly the same. The original four-fold 
symmetry enforced by the surface is broken because of the presence of the vacancy. For 
n = 3 and 4, the consequences of lifting the four- fold surface symmetry is that the t rimer 
is no longer an equilateral but a scalene triangle lying on the surface; the adsorption of 
the square structure now takes (instead of releasing) 0.67 eV. The square configuration is 
disfavored by more than 2 eV, compared to the tetrahedral configuration. In all cases, the 
sintering energy of Pd cluster next to an oxygen vacancy is lowered by 0.5 — 1.0 eV as 
compared to those on the stoichiometric surface. Particularly, Es is close to zero for the Pd 
dimer and t rimer. It can be concluded that the oxygen vacancy does not act as nucleation 
site but repels and dissociates Pd clusters if deposited directly on top of the vacancy. 

For the pure surface with an oxygen vacancy, strong electron localization is observed, 
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M=2; £^^=1.79eV 





M=0; £^^=2.36eV 



2.55 (d) 




M=0; =1.99eV 



FIG. 5. Adsorption geometries of Pd^ cluster for n = 1 — 4 nearby an oxygen vacancy. The 
corresponding magnetic moment M and adsorption energies E^d are given below the geometries. 
Selective bond lengths (in A) of Pd-Pd and Pd-0 bonds are shown with arrows. For the Pd 
monomer in (a), the electron localization sites are shown in golden color. The gas phase cluster 
structures are shown in the insets. 



and the spin-triplet states is 1.6 eV lower in energy than the spin-singlet. When attached 
with Pd clusters, the situation changes dramatically. Except the case of Pdi, where the 
electronic character of the surface vacancy is largely preserved because of the non-direct Pd 
adsorption, in all other cases the ferromagnetic coupling between the two localized electrons 
is attenuated. The localization on individual Ti ions is completely quenched because of the 
surface-cluster interaction. In fact, if the Pd monomer is moved directly above the vacancy 
site, similar attenuation of surface spin states is achieved. 



Except for the case of Pdi, the Bader charges in Tab. [TTl]show that the usual "electron- 
donor" role played by the Pd cluster is now reversed. Overall the surface vacancy is back 
donating about 0.4e to 0.5e to the cluster — a significant amount of charge to infiuence 
possible catalytic reactions taking place on top of the cluster. The charge fiow explains why 
the Pd cluster doesn't bind strongly to the underlying Ti atoms on the vacancy or stays near 
compared to the stoichiometric case. From Tab. [TT] we know that Pd cluster gets stabilized 



14 




E(eV) 



FIG. 6. DOS of the adsorbed Pd^ clusters (n = 1 — 4) close to an oxygen vacancy. The DOS 
of the pure defective surface is plotted in black dashed lines for comparison purpose. The spin- 
resolved DOS (up and down) with the adsorbates are shown as red solid and blue dash-dotted lines 
respectively. Insets: selected states with position on the DOS marked by arrows. The position of 
the oxygen vacancy is marked with "V" . 

by donating part of its electrons to the surface. In the case of the Pd monomer bound to 
the bridging site shown in Fig. [5| the charge transfer from the Pd atom to the surface is 
0.28e. However we found that when the monomer is moved closer to an oxygen vacancy, 
electrons flow in the opposite direction. For h^ sites, the charge accumulation on the Pd 
monomer can be as large as 0.50e, indicating that the localized electrons on the Ti ions have 
a strong tendency of reducing Pd, and eflFectively decreases the energy beneflt of the Pd 
cluster forming a bond with the underlying surface. Binding to atoms with high electron 
affinity can reduce such energy penalty, leading to the bridging site between two O5 atoms 
as the most stable adsorption site for the Pd monomer. 

The DOS of the adsorbed Pd^ clusters and selected gap states are shown in Fig. jH) It 
can be seen that the gap states have almost no overlap with nearby oxygen atoms. Unlike 
the case of adsorption on the stoichiometric surface, the back-donation of electrons from 
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the surface preserves the electronic structure of the cluster in gas phase. These additional 
charges also made the overall system favor the spin-singlet state. 

D. Pd^ on the surface with a Ti-interstitial atom 

Conventionally, the oxygen vacancy was believed to be largely responsible for the observed 
gap states (with some contributions from surface hydroxyld^^'^^'^^'^ . However, Wendt et al 
pointed out that under certain experimental conditions, to a large extent the surface chem- 
istry on reduced titania is associated with the Ti-interstitials existing in the near surface 
region.^ The additional Ti atom may be located at various depths.^ ^*^^*^^* Through intersti- 
tial exchange mechanism, these Ti-interstitial atoms can diffuse from the subsurface to the 
surface layer. To avoid potential artifacts associated with our computational setup where 
two bottom trilayers of Ti02 were fixed to their bulk positions, we only consider the case 
that a Ti-interstitial is located in the first trilayer, as shown in Fig. [7| 

The interstitial Ti atom reduces the surface with its four valence electrons, two of which 
are found to localize on the surface Ti^f atoms next to the interstitial atom. The other two 
localize in the subsurface layer. The magnetic coupling between these electrons is relatively 
weak, giving an energy difference of only 0.01 eV between the spin-triplet (M = 2 /llb) and 
the spin-quintet (M = A iib) states. The spin-triplet is slightly lower in energy. When a Pd 
cluster is attached to the surface, a preference over the spin-quintet was identified. In what 
follows, to facilitate direct comparisons of the electronic structures of clusters of different 
sizes on the equal footing, we show only the results of spin-quintet states. The magnetization 
and the qualitative picture of localization confirms the one obtained by Mulheran et al.^ who 
studied the diffusion of Ti interstitials with DFT+[/ and an atomistic potential method.^ 
However, we obtained splitted DOS with the HOMO placed 1.2 eV below the CBM, whereas 
in their paper all Ti3rf states share the same energy, whose distance to the CBM is less than 
0.5 eV. The HOMO and HOMO-1 are shown in the insets of Fig. [7| As one may see, though 
strongly localized in nature, the excess electrons are, in fact, shared by more than one Ti 
ions. 

To the best of our knowledge, the adsorption of Pd^ clusters (n = 1 — 4) on the titania 
(110) surface with a subsurface Ti-interstitial atom has not been investigated yet. The 
obtained geometries and adsorption energies are shown in Fig. [s] (a-d). In addition, the 
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FIG. 7. DOS of the Ti02 surfaces with a subsurface Ti-interstitial atom. The inset shows the 
geometry of the surface model. Ti ions with locahzed electrons are highlighted in golden color. The 
spin resolved DOS are shown as solid black and dashed red lines for spin up and down respectively. 
Selected gap states formed by localized electrons are shown in the insets with arrows marking their 
positions on the DOS. 



Bader charges and sintering energies are given in Tab. IV, The Pd monomer, again, strongly 
prefers to bind to the O5 atoms. The adsorption geometry is very similar to the stoichiometric 
result (see Fig. [2] (a)), though the adsorption energy Ead is 0.19 eV higher. Contrary to the 
case of oxygen vacancy, the presence of subsurface Ti-interstitial atoms may attract Pd 
monomers and act as nucleation centers for cluster growth. 

For the Pd dimer, similar to the other cases, multiple adsorption geometries with small 
energy differences were identified. The lowest energy configuration is shown in Fig. [s] (b). 
The corresponding spin-triplet state has a similar geometry but higher in energy by 0.58 
eV. Compared to other types of surface, the subsurface Ti-interstitial atom stabilizes the 
adsorption of the Pd dimer by 0.2 eV. The surface structure is strongly distorted owing to 
the adsorption of the Pd2 cluster. One Ti^f atom between the cluster and the interstitial 
atom is forced to move out of the surface layer and form a bond with an adjacent O5 atom. 
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FIG. 8. Adsorption geometries of Pd^ cluster for n = 1 — 4 with an Ti-interstitial. The magnetic 
moment M and adsorption energies Ead are given below the geometries. Selective bond lengths 
(in A) between Pd-Pd and Pd-0 are shown with arrows. 



This type of surface reconstruction is very similar to the minimum energy pathways obtained 
in the simulated diffusion of Ti interstitial atoms from subsurface to surface layer .^^^^ 

Pd3 and Pd4 show a similar tendency of binding to the protruding surface layer Tis/ 
atom (see Fig. [s] (c-d)). For Pda, we are unable to reach self-consistency for the spin-triplet 
state, whereas the spin-singlet state has the same energy as does the spin-quintet. In the 
case of Pd4, the situation is reversed, the calculation could not reach self-consistency for the 
spin-singlet. But the energy of the spin-triplet state is almost identical to that of the spin- 
quintet one. Due to the presence of the Ti-interstitial atom, the four-fold surface symmetry 
is effectively lifted. The binding geometries are no longer in-plane 2D but show a tendency 
of growing along the surface normal. 



The average Bader charge per Pd atom are shown in Tab.jTVj The charge transfer behavior 
shows an interesting sign change as n increases. The Pd monomer, unlike the case of oxygen 
vacancy, is an electron donor rather than an acceptor. The Pd trimer is slightly negatively 
charged. The Pd dimer and tetramer have unequally distributed charges over their atomic 
components, with an average charge transfer close to zero. Overall the charge transfer 
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n Geometry M (/is) Ead{(i\) Es{eY) AQ[Pd](e/atom) 



1 Monomer 


4 


1.72 




-0.20 


2 Dimer 


4 


2.55 


0.38 


-0.01 


3 Triangle 


4 


2.69 


0.86 


0.11 


4 Tetrahedron 


4 


2.15 


0.63 


0.00 



TABLE IV. Magnetic moment, adsorption energy Ead^ sintering energy Eg, and electron transfer 
AQ of Pd^ (n = 1 — 4) clusters adsorbed on Ti02 surface with a subsurface Ti interstitial atom. 

between the cluster and the surface is not significant. This may be related to the fact that 
both Pd and the Ti-interstitial atom are competing with each other to donate electrons to 
the surface. 

The spin-resolved DOS of the adsorbed Pd^ cluster and selected gap states are depicted in 
Fig. [9| Similar to the case of oxygen vacancy, there is no significant bonds forming between 
the cluster and the surface. Most gap states have similar character of Pd cluster states in 
gas phase, with its directions polarized by the underlying surface. For the case of the Pd 
monomer and dimer, the high magnetic moment is contributed largely by the surface. We 
obtain four Ti ions with one of its 3d states nearly filled by an electron. From the DOS, one 
can see that the diflFerence between the DOS of the up spin and down spin nearly coincides 
with the pure surface DOS. When n = 3 and 4, the cluster-surface interaction partially lift 
the electron localization — only three electrons are found localized on Ti ions. The remaining 
IjLiB contribution to the magnetic moment has no clear spatial origin. 



E. Cluster growth and sintering 

Clusters deposited on top of metal-oxides surfaces may absorb or lose a number of atoms 
and change their sizes at non-zero temperatures. Such sintering effects often lead to annihi- 
lation of small-sized clusters, which might be crucial for catalyzing desired reactions .^^^ 

The traditional theory on cluster sintering suggests that most cluster nucleation and 
growth processes involve the exchange of monomers between diflFerent clusters.'^ It is as- 
sumed that monomers are much more mobile than larger clusters. Using scanning tunneling 
microscopy (STM), Xu et al investigated the Pd cluster nucleation and growth on the 
titania surface experimentally.'® At low coverage (0.012 ML), they found that Pd clusters 
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FIG. 9. DOS of the defective Ti02 surfaces with subsurface Ti-interstitial atom. The inset shows 
the structure of the corresponding surface. The spin resolved DOS are shown as sohd red and 
dashed-dotted blue lines for spin up and down respectively. The DOS of the pure surface is plotted 
in black dashed line. The two gap states with localized electrons are shown in the insets with 
arrows marking their positions on the DOS. 



grow in 2D, with structures pseudomorphic to the substrate Ti02 surface. In their exper- 
iment, only the Pd dimer, t rimer and larger clusters were observed, and no existence of 
the Pd monomer was confirmed. On the other hand, Kaden et ai have measured the cat- 
alytic activity of CO oxidation of size-selected Pd cluster deposited on the titania surface. 
Clear diflFerences in the catalytic activities of Pdi and Pd2 were observed, suggesting the Pd 
monomer can exist on the surface at UHV condition. 

To quantify the mobility of the Pd monomer, we calculated PESs of a Pd atom adsorbed 
on the stoichiometric and the defective titania. On the stoichiometric surface, we confirm 
that the Pd monomer is extremely mobile along the O5 rows, forming a "two-way traffic" 



between two O5 rows with the h^ sites as energy minima. The PES is shown in Fig. 10 (b). 
The energy barrier between energy minima is only 0.05 eV. To cross the O5 row, the energy 
minimum path involves the Pd atom first moving to the site between two adjacent O5 atoms 
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then further moving to the other side. The highest barrier along the path is 0.3 eV, which 
is accessible at the room temperature. When the Pd monomer moves nearby an oxygen 
vacancy, the h^ sites are still the local energy minima, however the bridging sites between 
two O5 atoms one row away from the vacancy are more favored by the Pd monomer, as 



discussed in Sec. IIIC, The vacancy itself is very difficult to reach (energy barrier larger 



than 1 eV), though there is an energy minimum right on the vacancy once the barrier has 



been overcome (see Fig. 10 (c)). In the case that a Pd atom is close to a subsurface Ti- 
interstitial atom, the PES is more distorted by the presence of the additional Ti atom. 
We found a small spot for the Pd atom to bind to the underlying interstitial atom, with 
surrounding high energy barriers. The PES rapidly approaches the stoichiometric one when 



Pd moves away from the Ti-interstitial atom, as shown in the right half section of Fig. 10 (d). 

To simulate the thermodynamics of cluster growth and sintering effect, we used the 
PES of the Pd monomer on the three types of surfaces, and construct a 2D square lattice 
model describing the movement of the Pd monomers with the Monte Carlo method. The 
lattice sites are chosen to be a repeating two rows of h^ sites and one row of bridging sites 
between two O5 atoms. The Pd atom is allowed to jump to the nearest and the next nearest 
neighbors. The hopping rates are modified using the harmonic transition state theory,^'' in 
which the barriers are read from the saddle points between the adsorption sites on the PES. 
We assume the classical cluster sintering model, in which only the Pd monomer is mobile. 
The larger clusters grow or diminish by acquiring/losing the Pd monomers. Each time we 
randomly pick one Pd atom and move it according to the Metropolis algorithm. If one Pd 
atom collides into another Pd^ cluster, a Pd^+i cluster is formed with probability one. If 
the chosen atom is already in a cluster, an energy penalty of the corresponding sintering 
energy Es is added to determine whether that atom will leave the cluster. We simulate a 
16 X 16 super-lattice with the periodic boundary condition. Each cell of the lattice contains 
one Ti02 (4 x 2) unit, and is assumed to have a random type of surface. The PESs of 
the adjacent cells are assumed to be independent of their neighboring cells. To simulate 
typical experimental conditions, we chose 10% cells to have an oxygen vacancy and 5% cells 
with a Ti-interstitial atom. The remaining cells have the stoichiometric surface. The Pd 
coverage is fixed at 0.16 ML (^ 410 atoms) throughout. MC simulations were done over 
1000 initial configurations with randomly distributed Pd monomers. For each configuration, 
a multiple of 10^ MC runs at fixed temperatures were executed to ensure the system reaches 
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FIG. 10. Potential energy surfaces (PES) for the Pd monomer on the stoichiometric and defective 
surfaces, (a) Top view of the surface ceh (with a Ti-interstitial atom as an example) aligned with 
the PES graphs in (b-d). (b) PES on the stoichiometric surface, (c) PES on the surface next to 
an oxygen vacancy, (d) PES on the surface with a Ti-interstitial atom. (N.B., the discontinuity of 
(d) is due to technical reasons when the PES was generated.) The positions of the vacancy and 
the Ti-interstitial are marked by arrows. The color bar shows the corresponding energy: —Ead- 
For the dark red regions in (c) and (d), the energies are positive and exceed the value denoted by 
the color bar. 



its equilibrium distribution. Finally, our calculation is essentially zero-temperature based. 
Using molecular dynamics simulations, Miguel et al studied the stability of Pd clusters on 
the titania surface, and found that the surface-cluster interaction decreases appreciably for 
temperatures higher than 800 K^^^ Deposited clusters may adopt gas phase geometries 
(e.g. the tetramer may acquire a tetrahedral structure), and become metallic at higher 
temperatures. Such effects are not taken into account in our modeling. 

The density changes of the Pd dimer as a function of temperature on different types of 



surfaces are drawn in Fig. [TTJ The thermal stability of clusters of different sizes is generally 
determined by the relative strength of their sintering energies Es. At all temperatures, there 
is only a trace amount of the Pd monomers remaining at equilibrium, since the sintering 
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FIG. 11. Density C (in ML) of the Pd dimer (Pd2) as a function of temperature on different 
types of surfaces. The black sohd hne (Sto) shows the one on the stoichiometric surface. The blue 
dashed dotted line (Ov) gives the density on the surface with an oxygen vacancy, while the red 
dashed line (Int) corresponds to the surface with a Ti-interstitial atom. Three insets shows the 
density distribution of Pd^ cluster with respect to cluster size n at temperature 300 K, 700 K, and 
1000 K. The MC simulations were performed at 0.16 ML the Pd monomer coverage. 



energy Es of the Pd monomer is implied to be zero. At low temperatures, the most abundant 
cluster is the Pd dimer and, Pd^ clusters are generally immobile. The Pd dimer is very stable 
on the stoichiometric surface as well as close to a Ti-interstitial atom for temperatures below 
600 K. When the temperature reaches above 400 K, the possibility of finding a Pd dimer near 
an oxygen vacancy is practically zero. We found that this is also true for larger clusters, in 



line with the small sintering energies on the surface with an oxygen vacancy (see Tab. III). 
The sintering process in which the Pd dimers lose individual Pd atoms and exchanging 
these to grow into larger clusters begins to take effect at 600 K. The maximum of the size 



distribution moves to larger clusters, which can be told by comparing the insets of Fig. 11 



at different temperatures. The Eg of the Pd trimer is about twice that of the Pd dimer on 
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the stoichiometric surface and near a Ti-interstitial. This imphes that the trimer should be 
stable below 1200 K, which is partly confirmed by the size distribution diagram at 1000 K. 

IV. CONCLUSION AND ACKNOWLEDGEMENT 

We have performed a systematic analysis on the geometric and electronic structures of 
small Pd clusters on the stoichiometric and two defective Ti02 (110) surfaces. We find 
that unlike the noble metals, which prefer binding to the surface oxygen vacancies, the Pd 
clusters tend to adsorb on the stoichiometric surface or next to a subsurface Ti-interstitial 
atom. Owing to the surface-cluster interaction, Pd clusters have 2D in-plane structures 
when adsorbed on the stoichiometric surface. For defective surfaces, we demonstrate that 
DFT+C/ is capable of reproducing the experimentally observed gap-states and predicts that 
the localized electrons and the high-spin states of the surface may be quenched when a Pd 
cluster is adsorbed. We show that the Pd monomer is very mobile on the Ti02 surface and 
contributes greatly to the cluster sintering effect. 

The support form the ACS PRF grant 51052-DNI6 is acknowledged. Most calculations 
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